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1 Introduction 


Interest in SiC has resurfaced because of recent advances in material growth 
techniques. SiC has some material parameters that have an advantage over Si and 
GaAs devices for high frequency and high power applications. SiC has one of the 
largest reported bandgaps of common semiconductors (2.85 eV). Consequently, SiC 
has a breakdown field which is an order of magnitude higher than the breakdown 
field of Si or GaAs. The low field electron and hole mobilities of SiC are smaller 
than GaAs and Si mobilities, but the saturated electron velocity is higher. The 
thermal conductivity of SiC is approximately three times larger than the conduc- 
tivity of Si which means that larger amounts of heat can be dissipated successfully 
in SiC. These properties could lead to an excellent microwave and millimeter-wave 
power source. Some researchers have reported device quality SiC growth but the 
production at best remains very inconsistent. 1-4 Before appreciable resources are 
dedicated to material growth and characterization, it is important that a thorough 
study of the theoretical performance of SiC devices be carried out so that one may 
be able to better understand whether the dedication of resources is justifiable or 
not. A complete analysis of SiC devices will also permit a comparison between SiC 
devices and some of the existing device technologies. A preliminary comparison 
can be made by using some of the general figures of merit (for example, Johnson 
and Keys, see Ref. 1) but a more complete and accurate understanding can only 
come from individual device analyses. With this point in mind, the first concrete 
simulation results for SiC IMPATT diodes is reported here. The device operating 
characteristics and the typical power generating capabilites of SiC IMPATT diodes 



are studied at four different operating frequencies, namely, 10, 35, 60 and 94 GHz. 


2 Material Parameters 

Silicon carbide can grow in a large variety of polytypes with different structures 
and properties. The two most common polytypes are the a (hexagonal) and (3 
(cubic) polytypes. Some of the material properties important for IMPATT simu- 
lation have been measured for only one ploytype and not the other. Since there is 
no one standard method of growing SiC, most of the material parameters reported 
differ from one other because of the different quality of the material. It was de- 
cided to use the most promising experimental parameters that have been published 
and furthermore assumed that the parameters would not differ substantially from 
aSiC to /?SiC. There are still some parameters that have not been measured for 
any polytype. In this case a reasonable guess was made. These parameters will 
be measured in the future to further improve the accuracy of the simulation. An 
illustrative example of such a parameter is the hole saturated drift velocity which 
is important for double-drift IMPATT device analysis. Since there are no exper- 
imental results for this parameter, it was decided to use the same velocity as the 
one measured for electrons. Table 1 lists all the material parameters used and 
the references 4-7 where they were reported. An asterisk denotes that there are no 
published results and thus values were estimated. 

The field vs. velocity characteristic for the carriers, obtained from von Muench 
and Pettenpaul, 5 is expressed as: 



t ’(E) = v tat (l - e-*'*-) , 


( 1 ) 


where E is the electric field magnitude (in V/cm) and v Jot , /i are as given in Ta- 
ble 1. Because of the lack of information regarding the diffusion coefficients, the 
Einstein relationship was used to calculate the diffusion coefficients from the mo- 
bility measurements. A constant value for the diffusion coefficients was assumed. 
The ionization rates at room temperature for electrons and holes were obtained 
from Dmitriev et. al. 8 where different rates were reported depending on whether 
the applied electric field is oriented parallel to or perpendicular to the c axis of the 
crystal. It was decided to use only one case and the case where the electric field is 
perpendicular to the c axis was selected. The ionization rates used are as follows: 

0 P = 2.54 * loV* 11 * 10 '/*) (2) 


where a and 0 are the ionization rates for electrons and holes, respectively, and 
are given per centimeter. 

Table 1: Material parameters of silicon carbide (SiC) used for the simulation of double-drift 
IMPATT structures (*these are estimated values). 


Material Parameter 


Electron saturated velocity 
Hole saturated velocity 
Dielectric constant 
Electron mobility 
Hole mobility 
Electron diffusion constant 
Hole diffusion constant 
Thermal conductivity 


Value 



Reference 


von Muench and Pettenpaul 5 
* 

Patrick and Choyke 6 
Nishino et al. 4 
Nishino et al. 4 


Burgemeister et al. 7 






3 Device Design and Simulation 


Since it is expected that double-drift IMPATT structures can generate more 
power than their single-drift counterparts, the performance of double-drift IM- 
PATT structures was designed and computed. At each frequency the drift length 
wd was chosen as 


xjb d ** v * 


0.74 

2 / 


(4) 


where v is the particle saturated velocity (assuming equal values for electrons and 
holes) and / is the operating frequency. Since the same saturated velocity for holes 
and electrons is assumed, a symmetrical device with w n = w p will result. Next, 
it is assumed that wd comprises two thirds of the device length so that the total 
length w is 


w = -w D 


(5) 


The doping concentration was chosen so that at the desired operating current 
density the electric held at the two contacts would be approximately one third of 
the maximum electric field which occurs at the p-n junction. This design crite- 
ria does not have to be achieved rigorously but it should be made certain that 
the device is always in punch through. Operating current densities were chosen 
based on previous simulations of Si IMPATT devices at the same frequencies. All 
simulations were carried out at room temperature. 

Details of the device simulation algorithm were described by Bauhahn and 
Haddad . 9 The drift regions of the device are modeled by use of difference equation 
versions of the device equations where the carrier velocities, diffusion coefficients 
(assumed constant in this case), and ionization rates are field dependent. 



First, the dc solution was obtained and the doping level was varied until an 
appropriate electric field profile was acheived. For the ac solution, a sinusoidal RF 
voltage is impressed across the device so that the terminal voltage as a function of 

time is given by 

V(t) = V de + V RF sin(2irft) (6) 


Also, a particular dc bias current density./^ is imposed. Given these constraints, 
the program iterates until a periodic solution is obtained. At this point the device 
terminal current waveform is Fourier analyzed; the fundamental component of 
current together with the sinusoidal voltage component permits calculation of the 
device efficiency 


Prf 

P dc 


( 7 ) 


and generated power Prf at the fundamental frequency. The simulation results 
are normalized per unit area so that Prf is in W /cm 2 . 

Devices designed with the above-mentioned design rule were simulated at 10, 
35, 60 and 94 GHz. It is found that the device efficiencies are considerably lower 
than the efficiencies obtained from Si and GaAs IMP ATT devices at comparable 
frequencies. The lower efficiencies are not a result of some fundamental material 
limit, but they result because of a nonoptimum device design. Since the mobility 
of the holes is considerably lower than the mobility of the electoms, the field in 
the p region is not large enough to cause the holes to travel with the saturated 
velocity during the whole period of the drift time. Thus, at the end of the RF 
cycle, while all the electrons have been collected by the contact electrode, there 
are still holes (and thus charge) left in the device. This physical phenomena results 
in the degradation of the device efficiency. To correct this nonoptimum design, the 


p region of the double-drift device was made approximately half the value of the 
n region. Accordingly, the doping of the p region had to be changed so that the 
field at the contact was still approximately one third of the maximum field. In this 
way, the holes have to traverse a smaller distance than the electrons. Figure 1 
shows the schematic and electric field profile of such a device. Table 2 shows the 
doping and device dimensions calculated for non-symmetrical devices operating at 
four different frequencies. 

As in the case of the symmetrical devices, first a dc solution was obtained 
which was then used as the starting point for the large-signal ac solution. The 


+ + 




Figure 1: Doping profile and electric field profile for the simulated double-drift 

IMPATT structures. 





Table 2: Doping densities and device lengths at different frequencies for nonsymmetrical double- 
drift diodes. 


f (GHz) 

No (cm -3 ) 

tv. (H 

N a (cm- 3 ) 

W p ^m)_ 

10 

7.0xl0 19 

10.4 

1.3x10 18 

5.6 1 

35 

2.2x10“ 

3.10 

5.0x10“ 

1.55 

60 

5.0x10“ 

1.80 

9.0x10“ 

0.9 

94 

9.4x10“ 

1.18 

b-* 

X 

M 

o 

-4 

0.6 i 


results obtained from the ac simulation are presented in Tables 3 through 6. Each 
table lists results for a different operating frequency. The first five columns of 
the tables characterize the operating conditions of the IMPATT devices. For each 
case, a series of simulations are first presented where Vrf is varied but Jj c is fixed. 
Vrf is varied until a point is reached where the efficiency of the diode becomes 
maximum. Vrf is then held to this constant value and Jdc is now varied. Jdc is first 
decreased, then increased from the starting value for each case. The point in doing 
two sets of simulations is to try to optimize the device efficiency and find out the 
power generating capabilities of the diodes. However, this indeed does not assure 
an optimum design. In order to truly optimize the design, even more thorough 
calculations are required. However, the calculations presented do provide a typical 
and reasonably accurate picture of the device performance. 

Figures 2 through 5 show the large-signal results for the four devices operating 
at different fundamental frequencies. Table 7 gives a summary of the SiC IMPATT 
diodes when matched to a resistance of 1 0. From the simulated results presented 
in the table, it is noted that large power densities can be obtained at microwave 
and millimeter-wave frequencies with SiC double-drift IMPATTs. However, as is 
shown next, there are practical limitations that significantly decrease the power 
generating capabilities of these diodes, especially at millimeter- wave frequencies. 




Table 3: Operating properties of SiC nonsymmetrical double-drift IMPATT devices at 10 GHz. 


Vrf 
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Vic 

PRF 

Jdc 

Area 

Prf 


1 

1 1 

(Per- 
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1^91 




BSZEgB 

(kA / cm 2 ) 

_W) 

(kW-lfi) 



32.6 

10.2 

1836 

2.25 

1.20 

ig 



1.47 

32.9 

16.0 

1796 

3.60 

1.25 
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1.13 

33.1 

21.1 

1734 

4.60 

1.26 

1 

in 



33.2 

22.1 

1719 

4.65 

1.22 


0.435 


0.945 

33.3 

22.7 

1698 

4.73 

1.22 
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0.404 

1100 

0.713 

33.5 

21.8 

1645 

4.32 

1.20 
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0.273 

1200 

0.416 
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15.6 
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1.21 
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1000 


33.9 
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1000 

1 

33.9 
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1000 


33.8 

17.6 
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0.014 
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0.358 
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0.425 
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21.1 
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0.394 

1000 

1.24 

33.0 

22.8 

1695 

6.18 

1.60 
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3.42 

30.4 

20.4 
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5.04 
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3.90 
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6.17 
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8.71 

1000 

4.10 

28.3 

17.6 

1637 

20.5 

7.14 

5.010 

10.3 
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4.19 

27.2 

15.7 

1620 

21.0 

8.27 

5.550 

11.6 

1000 

3.93 

25.3 

12.3 

1594 

19.6 

10.0 

6.000 

11.8 

1000 

3.49 

23.2 

9.24 

1584 

17.5 

11.9 

6.320 

11.0 
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3.16 
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7.79 
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15.8 
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14.0 

14.1 
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8.52 

1000 
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20.4 
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1560 

11.0 

15.0 
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5.80 
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Table 4: Operating properties of SiC nonsymmetrical double-drift IMPATT devices at 35 GHz. 
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28.2 

22.0 
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425 
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31.3 
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30.7 
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55.4 

3.29 
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Table 5: Operating properties of SiC nonsymmetrical double-drift IMPATT devices at 60 GHz. 
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Table 6: Operating properties of SiC nonsymmetrical double-drift IMPATT devices at 94 GHz. 
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Figure 3: Large-signal results for a SiC double-drift IMPATT structure at f = 
35 GHz. (a) J = 3.46 kA/cm 2 , (b) Vrf = 425 V. 
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Table 7: Summary of properties of nonsymmetrical double-drift SiC IMPATT device performance 
when matched to a 1-fi resistance at different operating frequencies. 
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4 Power Limitations Of IMPATT Diodes 


There are two limitations that significantly reduce the power that can be gen- 
erated from IMPATT devices. First, there is the electronic limitation that comes 
from the restriction of matching the device impedance to the external circuit 
impedance and second, the thermal limitation that arises because of maximum safe 
operating temperature. For pulsed operation, if it is assumed that the duty cycle 
is so low that appreciable heating does not take place, the electronic impedance 
matching limitation would ultimately limit the maximum power generated by these 
devices. For cw operation, the limitations imposed by both considerations are cal- 
culated and the smaller value of the power is taken as the typical power expected. 

4.1 Electronic Limitation. The device admittance per unit area at the fun- 
damental frequency is given by Y = G + jB. Thus the device reactance can be 
calculated as 


X = 


B 

G* + B* 


( 8 ) 


and the resistance of the device is 








in Cl- cm 2 . To obtain the resistance in ohms, the preceding expression must be 
divided by the cross-sectional area (A) of the device. If it is assumed that practical 
circuit considerations limit the diode negative resistance to a minimum magnitude 
of 1 Cl, then the maximum allowed device area can be found such that (R/A )= 1 
and thus 

Q 

(G 2 + WjA = 1 ' (10) 

The impedance matching consideration sets a lower limit on the allowable mag- 
nitude of diode resistance and impedance. If the area were chosen larger, the 
negative resistance would be smaller than 1 Cl and Prf would be greater, but this 
will be difficult to match to any practical circuit. 

The device area and generated power for the case where the diode exhibits a 
negative resistance of 1 Cl are given in columns six and seven of Tables 3 through 
6. It is again noted that large amounts of power can be obtained. However, it 
should be remembered that these results are for pulsed operation where the heat 
generated is not a limiting factor. Table 8 gives a summary of the RF power 
obtainable from double-drift SiC devices at different frequencies. A comparison is 
also made with the power values obtained from Si and GaAs. As expected, SiC 
devices can generate significantly more pulsed power than Si or GaAs devices at 
all frequencies. 

Table 8: Comparison of peak electronic RF power (W) obtained at four different frequencies 
by matching the IMPATT device to a 1-0 load resistance (values for Si and GaAs diodes were 
obtained from Mains and Haddad*’ 10 ). 


Diode Type 

Frequency (GHz) Q 

10 35 60 94 | 

SiC uniform double drift 
GaAs uniform double drift 
Si uniform double drift 

11.8xl0 3 1.23x10 s 621 350 ) 

122 60.4 I 

24.3 I 






4.2 Thermal Limitation of IMPATT Diodes. A proper design requires 
that the device be able to dissipate the maximum allowed input power without 
exceeding a maximum allowed temperature. When the diode is operated in a 
continuous-wave (cw) mode, the thermal properties of the material become very 
important. In the case of SiC, since large amounts of power are theoretically possi- 
ble, thermal limitations can become very significant in determining the achieveable 
power densities during cw operation. 

In order to estimate the thermally limited power from the devices, the thermal 
resistance of the diode must be known. To calculate the thermal resistance of 
SiC IMPATT diodes an expression is used that has been used successfully for 
GaAs IMPATT diodes. 10 ’ 11 The only change made in the original expression is the 
substitution of the SiC thermal conductivity value for that of GaAs. For double- 
drift SiC IMPATT devices, the thermal resistance can then be given approximately 

by 

_ 2 _ 4l\T\ AljTj 

Rth “ + Rpk9 + 989*<P + 2322 ’ ( U) 

where d is the device diameter (in cm), Kh» is the thermal conductivity of the 
heat sink material (in W/cm-°K), Rpk g is the portion of the thermal resistance due 
to the package and bonding (in °K/W), /i is the buffer or substrate thickness (in 
cm), T\ is the average temperature of the buffer layer (in °K), / 2 is the active layer 
thickness between the junction and the buffer layer (in cm), and T 2 is the average 
temperature in the active layer (in °K). 

The first term is the spreading resistance that occurs at the interface between 
the diode and the heat sink. It is the result obtained by Kennedy 13 for the case 
of a uniform current heat flux density incident on a large cylindrical heat sink of 
thermal conductivity K^. The third term in the equation, the thermal resistance 
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of the buffer layer of length /j, can be derived from an expression for thermal 
resistance, assuming one-dimensional heat flow, given in Olson 13 

R tH = &x/KA , (12) 

where Ax is the length of the section, K is the thermal conductivity, and A is 
the cross-sectional area. Setting Ax = /j and equating these two terms yields 
the following expression for the thermal conductivity of SiC in the buffer layer (in 
watts per centimeter per degree Kelvin): 

K = 989/Tx . (13) 

The fourth term in the expression represents the thermal resistance of the active 
layer and is of the same form as the buffer layer term but is reduced by a factor 
of 0.4. This is because, in deriving the active layer term, uniform heat generation 
throughout the active region is assumed. Since R p k g is very dependent on the 
packaging and bonding techniques, it is taken to be zero in the thermal resistance 
calculations. For these calculations, the thermal conductivity of a diamond heat 
sink is taken as 11.7 W/cm-°K and that of a copper heat sink as 3.9 W/cm-°K. 

The use of multiple mesas or ring geometries will reduce the first term in the 
equation as explained in Mains and Haddad. 11 It is assumed here that an improve- 
ment of the spreading term by a factor of 0.55 is possible by utilizing geometries 
more efficient than the single round mesa structure. 

Once the thermal resistance expression is known, the diode junction tempera- 
ture T under large-signal cw conditions is given (in degrees Kelvin) by 


T = ZOO + R th (l - , 


(14) 



where R tfi is the thermal resistance, tj is the fractional RF generation efficiency, 
Vde is the diode operating voltage (in V), and Jdc is the diode current density (in 
A/cm 3 ). For this case it is assumed that /i = 1 /im and T\ = 725°K. /?, of course, 
depends on the device and was selected to be 750°K. 

From the earlier simulated results it is noted that as the current density is in- 
creased the power generated in the device also increases until it reaches a saturated 
value. However, higher currents also result in large amounts of power being dissi- 
pated inside the device which increase the junction temperature. To calculate the 
maximum amount of current that may be used, equations (10) and (13) are solved 
simultaneously and the temperature is set at T=773°K. It is assumed that this 
is the maximum allowable temperature for SiC devices. The equations are solved 
in terms of the diode diameter. This places a maximum limit on the area of the 
diode. 

Four different IMPATT cases are analyzed for thermal limitations. In the first 
case a copper heat sink is used with a single mesa structure. The second case 
considers a diamond heat sink with a single mesa said the third and fourth cases 
consider a ring geometry with copper and diamond heat sinks, respectively. 

5 CW Power Generation 

Two methods of calculating the maximum diode diameter and maximum power 
generated are now available. One is from impedance matching considerations in 
which a minimum circuit resistance of 1 Cl is assumed. The second method is from 
thermal considerations. To estimate the maximum achievable power in practice, 
the procedure is to calculate the device diameter using both methods and then to 
choose the smaller of these two values. This procedure has been carried out for 


the devices under consideration and the maximum achievable cw power is given in 
Table 9. Table 10 shows a comparison of SiC, Si and GaAs diodes operating in the 
cw mode at different frequencies. The performance of SiC at lower frequencies is 
superior to GaAs and Si IMPATT performance. However, at millimeter-wave fre- 
quencies, Si devices can produce similar amounts of cw power as the SiC IMPATT 
devices. 


Table 9: Summary of continuous wave (cw) RF power (W) for SiC double-drift devices at four 
different frequencies for mesa and ring structure geometries and different heat sink materials. 


Diode type 

Frequency (GHz) 

10 35 60 94 

Single mesa, copper heat sink 
Single mesa, diamond heat sink 
Ring structure, copper heat sink 
Ring structure, diamond heat sink 

5.14 2.61 1.74 0.55 

28.1 11.0 6.03 2.57 

11.9 5.82 2.76 1.58 

56.20 22.0 11.0 5.23 


Table 10: Comparison of cw RF power (W) for various double-drift IMPATT devices with two 
different heat sink materials and two different geometries (values for Si and GaAs diodes are from 
Mains and Haddad 9 ' 10 ). 


Diode 

Material 

u 

Single Mesa 

Ring Structure 

Copper Sink Diamond Sink 

Copper Sink Diamond Sink 

Si 

GaAs 

SiC 

10 

0.967 3.93 

5.14 12.2 

2.32 5.91 

30.3 56.2 

Si 

GaAs 

SiC 

35 

2.61 11.0 

5.82 22.0 

Si 

GaAs 

SiC 

60 

■m ■■ 

0.63 1.8 

2.76 11.0 

GaAs 

SiC 

94 


0.099 0.181 

1.58 5.23 




















There are two points that should be noted. First, in the thermal resistance 
calculations it is assumed that the device is operating at 773°K while the results 
of the room temperature simulations are used. This was done because there is 
not enough information to correctly predict the material parameters at 773°K. 
Since, in general, device performance degrades with increasing temperature, the 
cw results are perhaps a bit optimistic. Second, the electric field in these structures 
gets extremely high and there might be some possibility of tunneling. This was 
not incorporated into the simulation. Again, this could result in different device 
operating conditions. 

6 Conclusions 

SiC double-drift IMPATT diodes were simulated at 10, 35, 60 and 94 GHz using 
available material parameters. In the pulsed mode, the SiC IMPATT performance 
was found to be far superior to the Si and GaAs IMPATT performances. In the 
cw mode, the performance of SiC is comparable to Si double-drift devices but 
superior to GaAs devices especially at lower frequencies. Therefore, significantly 
improved power levels may be expected from SiC double-drift IMPATTs used in 
the pulsed mode, but thermal considerations severely limit the power available for 
cw operation at millimeter-wave frequencies. 
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Abstracts 


Microwave and millimeter-wave power generation in silicon carbide 
avalanche devices 
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University of Michigan. Ann Arbor, Michigan 48109-2122 

(Received 3 December 1987; accepted for publication 15 March 1988) 

Silicon carbide (SiC), due to its thermal and electronic properties, has long been considered an 
excellent device material for microwave and millimeter-wave power generation. Numerical 
simulations were performed to study the typical power generating capabilities of SiC impact 
avalanche transit-time (IMPATT) diodes utilizing the recent experimental data available. 
Operating characteristics of double-drift IMPATT devices at 10, 35, 60 and 94 GHz are 
compared. Both pulsed mode and continuous-wave (cw) mode operation are studied. Finally, 
a comparison among SiC, Si, and GaAs double-drift IMPATT devices is made at various 
frequencies. It is shown that, for the pulsed mode of operation, SiC double-drift IMPATT 
devices can produce significantly higher powers than Si and GaAs devices at comparable 
frequencies. In the cw mode of operation, SiC devices can produce significantly more power 
than GaAs devices at all frequencies. However, a comparison at 94 GHz indicates that SiC 
IMPATT diodes in the cw mode of operation produce power levels comparable to Si IMPATT 
devices. At lower frequencies the performance of SiC diodes operating in the cw mode is 
expected to be better than the performance of Si devices due to the better thermal conductivity 
of SiC. 
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